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ABSTRACT 

Using the results from a previously developed Lya/continuum produc- 
tion/transmission and dust enrichment model for Lyman Alpha Emitters (LAEs), 
based on cosniological SPH simulations, we assess the detectability of their dust- 
reprocessed sub-millimeter (submm) radiation. As supernovae (rather than evolved 
stars) control dust formation and destruction processes, LAEs are relatively dust- 
poor with respect to local galaxies: they have low dust-to-gas ratios (0.05 times 
the dust-to-gas ratio of the Milky Way) in spite of their relatively high metallicity, 
Z « 0.1 — 0.5^0. Using the derived escape fraction of ultraviolet (UV) continuum 
photons we compute the UV luminosity absorbed by dust and re-emitted in the far 
infrared. The LAE submm fluxes correlate with their Lya luminosity: about (3%, 1%) 
a.t z = (5.7, 6.6) of the LAEs in our simulated sample (those with with LogLa > 43.1) 
would have fluxes at 850 /im (the optimal band for detection) in excess of 0.12 mJy 
and will be therefore detectable at 5a with ALMA with an integration time of only 1 
hour. Such detections would open a new window on the physical conditions prevailing 
in these most distant galaxies. 
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1 INTRODUCTION 

Lyman Alpha Emitters (LAEs) are galaxies identified by 
means of their strong Lya emission. Unambiguous charac- 
teristics of this line including the strength, width, and blue- 
ward cut-olT have enabled recent detections of LAEs out to 
z ~ 7.7 (Hibon et al. 2009), adding to the already existing 
data at z ~ 7 (lye et al. 2006), z ~ 6.5 (Kashikawa et al. 
2006), and lower redshifts. Since LAEs are among the most 
distant objects known, in addition to being superb probes of 
reionization, they serve as laboratories to study early galaxy 
evolution and clarify fundamental but poorly understood as- 
pects of high-2 galaxies such as their dust content and ex- 
tinction law. 

Although a number of studies have been conducted, the 
extent and sources of the dust enrichment of LAEs remains 
a much debated issue. Using cosmological SPH simulations, 
Dayal et al. (2009a) and Nagamine et al. (2008) have shown 
that at z ~ 5.7, the color excess of LAEs, E{B — U) ~ 0.15. 
These theoretical estimates are not inconsistent with recent 
experimental determinations: by fitting the SEDs of 3 LAEs 
at z = 5.7 by Lai et al. (2007) have inferred E{B - V) < 
0.225-0.425; in a sample of 12 LAEs a.t z = 4.5, Finkelstein 
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et al. (2009a) have found A1200 = 0.5 — 4.5; finally, Pirzkal 
et al. (2007) have found Av = 0.05 — 0.6 for 3 galaxies at 
2 = 4 — 5.76. The latter two values, when translated into 
the color excess (using a supernova dust extinction curve, 
e.g. Bianchi & Schneider 2007) are found to be E{B -V) = 
0.035 - 0.316 and 0.025 - 0.3. These values therefore result 
in a coherent picture. However, Gronwall et al. (2007) and 
Ouchi et al. (2008) find much lower values of the color excess, 
E{B - U) s; 0.05 for z = 3.1 LAEs. 

The main uncertainty in using the SEDs to infer the 
dust content lies in the fact that the observed spectra de- 
pend on: (a) stellar population properties such as the age, 
metallicity, initial mass function (IMF) and (b) the amount 
and distribution of dust in the interstellar medium (ISM). 
These two effects are not easily disentangled from each other 
and this might be the source of the disparity in the dust ex- 
tinction inferred by different works. These difficulties can 
be overcome if one can exploit the fact that dust absorbs 
ultraviolet (UV) radiation and re-emmits it at far-infrared 
(FIR) wavelengths; thus, if LAEs could be observed in the 
sum- millimeter (submm; FIR in the LAE's rest frame), a 
new window on the study of LAE properties would open, 
aUowing to precisely determine their dust content and a far 
better determination of their stellar populations and star 
formation history. 
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Clumped dust has also been invoked by many authors 
(Dawson et al. 2004; Dayal et al. 2008, 2009a, 2009b; Finkel- 
stein et al. 2008, 2009a, 2009b; Kobayashi et al. 2007, 2009) 
to enhance the observed Lya equivalent width (EW) ; the es- 
cape fraction of Lya photons increases as compared to that 
of the continuum photons because of the so called 'Neufeld 
effect' (Neufeld 1991), if dust in the ISM is clumped. If how- 
ever, LAEs are found to contain negligible amounts of dust 
from FIR observations, scenarios which are considered ex- 
otic at present, including a top heavy IMF, PopIII stars and 
LAEs being extremely young objects (ages ^ 10 Myr) would 
have to be considered quite seriously. 

The main aim of this paper is to assess the detectability 
of LAEs in the submm band. We plan to accomplish this 
task by building upon a previously developed detailed model 
based on cosmological numerical simulations. Such model 
(Dayal et al 2009b) successfully predicts the Lya and UV 
Luminosity Function (LF) of the LAE galaxy population; 
in addition, it provides a wealth of information about the 
physical properties of these objects. 



2 METHOD 

2.1 Basic simulations 

We start the calculation from a previously developed model 
of LAEs, the details of which can be found in Dayal et 
al. (2009b). In brief, the cosmological SPH simulatiouii^ 
of (75/i~^Mpc)^ (comoving) have been carried out using 
a TreePM-SPH code GADGET-2 (Springel 2005) with the 
implementation of chemodynamics as described in Torna- 
tore et al. (2007). Galaxies are recognized as gravitation- 
ally bound groups of star particles by running a standard 
friends-of-friends (EOF) algorithm, decomposing each EOF 
group into a set of disjoint substructures and identifying 
these. After performing a gravitational unbinding procedure, 
only sub-halos with at least 20 bound particles are consid- 
ered to be genuine structures, Saro et al. 2006). For each 
"bona-fide" galaxy, we compute the mass-weighted age, the 
total halo/stellar/gas mass, the star formation rate (SER), 
the mass weighted gas/stellar metallicity, the mass- weighted 
gas temperature and the half mass radius of the dark matter 
halo. 

We have coupled the above simulations with a Lya/ 
continuum production/ transmission model as in Dayal et 
al. (2009a). For each galaxy, the total Lya and continuum 
luminosity are calculated as a sum of the contributions from 
(a) stellar processes and (b) cooling of coUisionally excited 
H I in the ISM (Dayal et al. 2009b) . In the post- processing, 
assuming Type II supernovae (SNII) to be the primary dust 
factories for z > 5.7, we include a model to calculate the to- 
tal amount of dust in each galaxy depending on its intrinsic 
properties such as the SFR, age, metallicity and gas mass. 
This dust mass is then used to calculate the associated op- 
tical depth, which is translated into the escape fraction for 



^ The adopted cosmological model for the simulation corresponds 
to the ACDM Universe with = 0.26, = 0.74, Ob = 0.0413, 
rts = 0.95, Ho = 73 km s~^ Mpc~^ and erg = 0.8, consistent with 
the 5-year analysis of the WMAP data (Komatsu et al. 2009). 



continuum photons, /c, assuming a slab-like dust distribu- 
tion, to obtain the UV luminosity function (LF). The escape 
fraction of Lya photons is then assumed to scale with fc to 
obtain the observed Lya LF. 

Galaxies are identified as LAEs based on the current 
operational definition, i.e. observed Lya luminosity. La ^ 
10^'^'^ ergs~^ and observed equivalent width, EW ^ 20 A. 
In total, (39000,18000) galaxies have been identified in the 
two relevant simulation outputs ai z — (5.7, 6.6) out of which 
(1045, 502) have been selected as LAEs according to the 
above criterion. 



2.2 UV escape fraction and FIR luminosity 

We estimate the total FIR luminosity for all the galaxies 
identified as LAEs at the redshifts of interest, i.e. z 
(5.7,6.6). For this calculation, we assume that the dust 
is predominantly heated by UV radiation from stars with 
wavelengths 912-4000 A, Buat & Xu (1996); in this paper, 
the term UV is used for non-ionizing continuum in this wave- 
length range. We start by calculating the optical depth to 
dust, Tuv, seen by UV photons as ruv = 3E(i[4as]~^, where 
Ed = Mdust / i'r^r'^) is the dust surface mass density, a and s 
are the grain radius and material density. The dust surface 
mass density is estimated by assuming that the dust mass, 
Mdust is spread on a scale, — (0.6, 1.0)re at z ^ (5.7, 6.6), 
where re is the stellar distribution scale inferred using the 
results of Bolton et al. (2008) . Complete details of this cal- 
culation can be found in Dayal et al. (2009b). 

Assuming the dust and stars to be homogeneously dis- 
tributed, we use this optical depth to calculate the UV lu- 
minosity, Luv, escaping from the galaxy as Luv = L'^vfc, 
where Lliy is the intrinsic UV luminosity, calculated using 
the population synthesis code STARBURST99 (Leitherer 1999) 
using the appropriate IMF, star formation rate (M*), metal- 
licity [Z) and age (t,) for each LAE, as obtained from the 
simulation. Further, fc, the escape fraction of the UV con- 
tinuum is calculated assuming a slab-like dust distribution, 
such that 



1 - e~ 



Tuv 



(1) 



Hence, /c depends on the intrinsic properties of the galaxies, 
i.e, the dust mass and dust distribution scale. The average 
value of {fc) = (0.23,0.38) at z ~ (5.7,6.6) respectively, 
Dayal et al. 2009b). 

Because of the large cross section of dust against UV 
light and the intense UV field in a star forming galaxy, we 
can assume that the emitted FIR luminosity, Lpjn, is equal 
to the UV luminosity which is absorbed and heats up the 
dust, such that 



LpiR = L%Y ~ Luv = (1 — fc)LijY- 



(2) 



In order to examine the detectability of dust emission, 
the observed fiux is predicted to be 



(1 + z)L^(x+z) 



(3) 



where dL is the luminosity distance (Carroll, Press & Turner 
1992), and is the monochromatic luminosity at the cho- 
sen frequency of observation. 
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If dust grains emit thermally with a single dust temper- 
ature, Tdust, Lv can be written as 



Lv — 'iTvMdu3tHuB^(Tdust), 



(4) 



where k^, is the mass absorption coefficient, B^{Tdust) is the 
Planck function with frequency f and temperature Tdust- 

If we assume a power law form for Ki, such that Ki, = 
K^g{u/uo)'^ , the total FIR luminosity can be expressed as 



-f3 fkTdustY+'' (2h 



dx, (5) 



where, k is the Boltzmann constant, h is the Planck con- 
stant, c is the light speed and x = hv{kTdust)~^ ■ We use 
vq — 3.00 X 10^^ Hz (corresponding to a wavelength of 100 
^m), which is an arbitrary frequency in FIR used for nor- 
malization. 

Since we assume all the dust grains to be spherical with 
a single size a — 0.05/im (appropriate for a SN-produced 
dust, see Todini & Ferrara 2001; Nozawa 2003, 2007) and 
material density s, the mass absorption coefficient can 
be written as n^g = 3Qv[4as]~^ , where is the optical 
absorption cross section normalized to the geometrical cross 
section (na^). For graphite/carbonaceous grains, Q^a"^ — 
1.57 X 10^ cm~^ at a wavelength of 100 fiva, s = 2.25 g 
cm~^ and /3 = 2 (Draine & Lee 1984). This results in k^, — 
Even though in this work, we assume 
the dust grains to be graphites, the flux is insensitive to the 
specific grain material for a given value of Lfir since a small 
Kv is compensated by a high Tdust and vice versa. Then, 
equation (O can be simplified to the following expression 



: 6.73 



Mdust/Mo 



1/6 



K. 



(6) 



This temperature is put in eq. Q to obtain Lv and hence 
the dust emission flux detectable in the submm bands in the 
observer's frame. 



3 RESULTS 

3.1 Dust abundance 

The relation between dust-to-gas ratio [D) and gas metal- 
licity is useful to discuss the dust enrichment scenario in the 
low-metallicity phase as shown by Lisenfeld & Ferrara (1998) 
for nearby blue compact dwarf galaxies. We calculate this 
ratio as _D = Mdust/Mg, where the dust mass, Mdust is cal- 
culated as shown in Dayal et al. (2009b). The gas mass, Mg 
and the mass- weighted gas metallicity, Zg, are both obtained 
from the simulation as mentioned before. We normalize the 
dust-to-gas ratio to the one measured in the solar neighbor- 
hood, Dq = 1/150 = 6 X 10"^ (Hirashita & Ferrara 2002). 
The relation between D and Zg has been investigated by sev- 
eral authors. In particular, Lisenfeld & Ferrara (1998) have 
shown that for nearby dwarf galaxies the following relation 



^ The corresponding values for silicates are Q^a ^ = 1.45 X 
lO'^ cm~^, s = 3.3 g cm~^ and /3 = 2, which results in fti/ = 
32.9(i^/uo)/^ cm2 g-i. 
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Figure 1. Relation between dust-to-gas ratio and gas metallicity 
for LAEs at z = 5.7 (circles), 6.6 (triangles). The dashed line 
corresponds to Zg oc D^-^-^ (Lisenfeld &; Ferrara 1998). 



holds: Zg oc D'', with p = 0.52 ± 0.25. In principle, there 
is no reason to expect that such low-redshift determination 
applies equally well to LAEs, due to the very different phys- 
ical conditions and shorter evolutionary timescales allowed 
by the Hubble time at 2: ~ 6. Using the above relation for 
a gas metallicity value equal to 0.2Zq, appropriate for our 
calculation of LAEs, yields D = 0.003 — O.I2D0; from our 
result we get a mean value in reasonable agreement with 
this expectation, D = 0.001 - 0.045Dq for Zg = 0.2Zq. 

However, although we confirm a correlation between 
dust abundance and metallicity, the D — Zg relation we 
find, shown in Fig. 1, shows considerable deviations from 
the above power-law. First, the slope of the relation is con- 
siderably steeper (i.e. p > 1); second, in spite of the rela- 
tively high metallicity of LAES (about 0.1 — 0.5^©), the dust 
abundance is relatively depressed with essentially all objects 
having D < 0.0573q. Taken together, these two points im- 
ply that LAEs are relatively poorly efficient dust produc- 
ers with respect to local galaxies. The physical explanation 
of this fact is straightforward. Because of their young ages 
(< 500 Myr), LAEs can only have their dust produced by 
SNII rather than by evolved stars, which dominate the dust 
production in local galaxies. However, since dust is both 
produced and destroyed by SN in LAEs, D saturates at a 
relatively low value; in local galaxies, instead, SN shocks 
cannot counteract the rate of dust production by evolved 
stars and D grows to larger values. 

Finally, Fig. 1 shows the presence of a considerable num- 
ber of outliers to the left of the power-law relation, i.e. the 
low D-high Zg region, at both redshifts. Since both the dust 
and gas mass scale with the halo mass (Dayal et al. 2009b) , 
this physically means that even some galaxies with a very 
small halo mass retain a large metal fraction in the ISM 
gas. This can possibly be explained by the fact that though 
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Figure 2. Correlations between Lya luminosity and submm flux for z = 5.7 (left) and 6.6 (right). The horizontal solid, dotted and 
dashed lines show the 5cr detection limit of ALMA with 1 hour integration for the bands corresponding to 100 Ghz (3mm), 220 Ghz (1.4 
mm) and 353 GHz (850/im) respectively. Symbols (squares, triangles, circles) in both panels show the results obtained from this work 
for the same bands (100, 220, 353 GHz). 



galaxies with a small halo mass have very feeble star for- 
mation rates and hence produce less amount of metals as 
compared to larger galaxies, the mechanical feedback that 
would eject their ISM metal content is comparatively more 
depressed with respect to larger galaxies. 



3.2 Prediction for sub-millimeter fluxes 

ALMA could be a strong tool to directly detect the dust 
emission from a large number of LAEs. Then, the question 
of how to identify the sources of submm emission arises. We 
find that the observed Lya luminosity at both z ~ (5.7, 6.6) 
(as calculated in Dayal et al. 2009b) and the fluxes at 3mm, 
1.4 mm and 850/xm of LAEs are correlated, as shown in Fig. 
[21 in general, the galaxies with the largest observed Lya 
luminosities also show the largest value of the FIR fluxes. 
The reason for this can be explained as follows: the observed 
Lya luminosity. La is related to the intrinsic value, L^"' as 

La ~ La faTa, (7) 

where fa is the fraction of Lya luminosity which emerges 
from the galaxy, undamped by dust and Ta is the fraction 
of Lya luminosity transmitted through the IGM. Now, the 
largest galaxies have the largest star formation rates (Dayal 
et al. 2009a) and hence the largest values of L^"* since this 
scales with the star formation rate. As shown in Dayal et 
al. (2009b), their Ta is also the largest. Further, the largest 
galaxies are also the most dust-enriched by SNII since the 
SNII rate scales with the star formation rate. However, due 
to large radii of dust distribution (which scales with the 
SFR), the largest galaxies do not have the smallest fa as 
expected. This implies that the largest galaxies show both 
the largest La and fi,, hence the trend seen. Fig. [2] shows 
that the 850/xm band is the optimum one to look for submm 
emission; (32, 5) LAEs at z = (5.7, 6.6) in our simulation 
volume have an 850/^m flux larger than the 5a 1 hour inte- 
gration limit of ALMA which is 0.12 mjy (Morita & Hold- 



away 2005). This corresponds to about (3%, 1%) of the LAEs 
in our simulated samples at these redshifts. 

In Fig. [S] we show the FIR luminosity and submm flux 
luminosity functions for our LAE sample at 850/im since this 
seems to be the most optimum band to look for FIR emis- 
sion, as mentioned above. We see that selecting LAEs with 
La S5 lO^^'^ergs"^, can efficiently detect galaxies whose 
submm fluxes are bright enough to be detected by ALMA 
with an integration time of only 1 hour. A longer integra- 
tion time drastically increases the number of LAEs detected 
by ALMA because of the steep rise of distribution function 
toward faint submm fluxes as shown in Fig. [3] 

Thus, we show that we can simply point ALMA towards 
bright LAEs that have already been identifled at high-z to 
detect LAEs in the sub-mm band at 5a with an integration 
time of only 1 hour. 



4 SUMMARY AND CONCLUSIONS 

This work is based on the results obtained in Dayal et al. 
(2009b), where, we coupled a cosmological SPH simulation 
with a Lya/continuum production/transmission model and 
a dust model. This enabled us to calculate the observed 
Lya/continuum luminosities and the dust mass evolution for 
each galaxy. We thus obtained the escape fraction of contin- 
uum photons for each galaxy based on its intrinsic properties 
including the star formation rate, age, metallicity, IMF and 
gas mass. The only free parameter, the escape fraction of 
Lya photons was obtained by comparing the model results 
to the observations at z ~ (5.7, 6.6). 

Though we find that the dust to gas ratio values as 
a function of gas metallicity obtained from our model are 
broadly consistent with those predicted by Lisenfeld & Fer- 
rara (1998) at Z = 0.2Zq, the slopes predicted by the two 
works are different. The reason for this disagreement arises 
because we assume SNII to be the only sources of dust. 
While this leads to a progressive metal enrichment of the 
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Figure 3. The distribution function of the observed 850^tm flux 
of LAEs at 2 = 5.7 (sohd line) and 6.6 (dotted line). The vertical 
dashed line shows the 5cr 1 hour integration detection limit of 
ALMA (= 0.12 mjy) and the shaded regions show the poissonian 



ISM and hence a high value of the gas metallicity, since 
SNII also destroy the dust that they produce, the dust to 
gas ratio remains very small (< 0.05 of that in the solar 
neighborhood). On the other hand, in dwarf galaxies, both 
SNII and evolved stars produce dust; the supernovae shocks, 
honrever, are not capable of destroying the huge amounts of 
dust produced by evolved stars and this leads to much larger 
values of the dust to gas ratio. 

For each galaxy identified as a LAE as these redshifts, 
we assume that the UV luminosity between 912-4000 A that 
is absorbed inside the galaxy heats up the dust, leading to 
re-radiation in the submm band in the observer's frame. This 
allows us to calculate the emitted flux at 3000, 1363 and 850 
/xm. 

We then present a very efficient strategy to look for 
dust emission from LAEs. We find that the observed Lya 
luminosity and the FIR fluxes are correlated for LAEs in 
all the three ALMA bands (3000, 1363 and 850 ^im); the 
LAEs with the largest observed Lya luminosity also show 
the largest fluxes. This means that pointing surveys of the 
brightest LAEs already identifled would yield large fluxes 
in ALMA. We find that the 850/im band seems to be the 
optimum one for detecting the FIR emission from high-z 
LAEs. Although a few tens of objects would be detectable 
with a 1 hour integration, the number density of objects 
detectable would rise steeply with the integration time. 

Finkelstein et al. (2009b) have used the measurements 
of the UV spectral slopes to derive far-infrared flux pre- 
dictions for a sample of 23 LAEs (Finkelstein et al. 2QQ9a; 
Pirzkal et al. 2007) at 4. It is interesting to compare 
to what extent the results obtained using their approach 
match with those from our model. We first compare the 
color excess obtained from these works where we find that 
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the average value, E{B — V) ~ 0.15 obtained from our 
work (Dayal et al. 2009a; 2009b) is very consistent with 
E[B -V)= 0.035 - 0.316 and 0.025 - 0.3, found by Finkel- 
stein et al. (2009a) and Pirzkal et al. (2007) respectively. 
However, while Finkelstein et al. (2Q09b) find that 39 ±22% 
of their LAEs would be visible with ALMA with an integra- 
tion time of 4 hours, we find only ~ (3%, 1%) of the LAEs 
in our simulation would be visible with a 5 — ct, 1 hour in- 
tegration limit at 2: ~ (5.7, 6.6). Although these results may 
seem discrepant at the first glance, they can be explained 
quite easily considering the redshifts of the LAEs modelled 
by Finkelstein et al. (2009b). In their sample, 20 of the LAEs 
are at 4 z ^ 5, of which, 9 (i.e. 45%) would be visible with 
ALMA with a 4 hour integration time. However, none of the 
3 LAEs with z > 5.2 would be detectable in the submm. 
This is highly consistent with our results, where only 3% 
of the LAEs would be visible with ALMA at z ~ 5.7 for a 
5 — (T 1 hour integration time; hence, in a sample of 3 LAEs, 
none would be visible. This is also a confirmation of our dust 
model, where SNII are the primary dust factories at z ^ 5.7, 
with evolved stars dominating at lower redshifts. Since SNII 
both produce and destroy dust, the dust enrichment levels 
at high redshifts are very small as compared to that at lower 
redshifts (2 5.7). Hence, the detectability of LAEs in the 
submm decreases with increasing redshift. 

The main uncertainty in this work lies in the fact that 
we assume SNII to be the primary sources of dust produc- 
tion/destruction at the redshifts considered. While this is 
justified for ages ^ 1 Gyr, we can not completely rule out 
the contribution from evolved stars as shown by Valiante et 
al. (2009). Further, the dust yield per SNII and the dust de- 
struction efficiency due to supernova shocks are only known 
well enough to within a factor of a few. Also, since the scale 
length of dust distribution is not well known at these red- 
shifts, we have to resort to the data to guide us. We plan to 
deal with these uncertainties in future works. 
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